Detection of brown dwarf-like objects in the core of NGC 3603 
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ABSTRACT 

We used near-infrared data obtained with the Wide Field Camera 3 (WFC3) on the Hubble 
Space Telescope to identify objects having the colors of brown dwarfs (BDs) in the field of 
the massive galactic cluster NGC 3603. These are identified through a combination of narrow 
and medium band filters which span the J and H bands and are particularly sensitive to the 
presence of the l.S-l.S/um H2O molecular band unique to BDs. We provide a calibration of 
the relationship between effective temperature and color for both field stars and for BDs. This 
photometric method provides effective temperatures for BDs to an accuracy of ±350K relative to 
spectroscopic techniques. This accuracy is shown to be not significantly affected by either stellar 
surface gravity or uncertainties in the interstellar extinction. We identify nine objects having 
effective temperature between f 700 and 2200 K, typical of BDs, observed J-band magnitudes 
in the range f9.5-21.5, and that are strongly clustered towards the luminous core of NGC 3603. 
However, if these are located at the distance of the cluster, they are far too luminous to be normal 
BDs. We argue that it is unlikely that these objects are either artifacts of our dataset, normal 
field BDs/M-type giants or extra-galactic contaminants and, therefore, might represent a new 
class of stars having the effective temperatures of BDs but with luminosities of more massive 
stars. We explore the interesting scenario in which these objects would be normal stars that have 
recently tidally ingested a Hot Jupiter, the remnants of which are providing a short-lived extended 
photosphere to the central star. In this case, we would expect them to show the signature of fast 
rotation. 

Subject headings: Stars: brown dwarfs - Stars: pre-main-sequence - Galaxy: open clusters and associa- 
tions: individual (NCG 3603) - Instrumentation: photometers (HST-WFC3) - Techniques: photometric 
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1. Introduction 

Understanding the origin of brown dwarfs 
(BDs) is an important component of the theory of 
star formation and evo lution that still r emain s 
somewhat mysterious ( Whitworth et all 20071) . 



Moreover, recent observations (Gould et al 



20091 ) 



suggest that old substellar field objects may be 
more common than previously assumed, opening 
again the debate on whether BDs are ubiquitous 
and occur in sufficient numbers to account for an 
appr eciable amoun t of dark matter in the Galaxy 
(e.g.. lTinnevlll999t) . 

Canonical approaches to these research fields 
include surveys for BDs in the field, studies of 
the sub-stellar initial mass function (IMF) in 
young clusters and associations under different 
star forming conditions, measurements of the disk 
frequency in the very low-mass regime and in- 
vestigations of the mass accretion process near 
and below the substella r boundary (see, e.g., 
Javawardhana et al.ll2003l) . However, these overall 
approaches ultimately rely on the spectral classi- 
fication of low-mass stars and BDs and, hence, on 
optical/near-infrared (IR) spectroscopy ( see, e.g.. 



McLean et all 120031 iGeballe et all 120021) . which 



has proved impossible for very faint distant BDs 
with current telescopes. Thus, the studies con- 
ducted so far are limited to galactic clusters and 
star forming regions within a few kpc from the 
Sun, while the substellar population of extreme 
environments, such as the outer/center of the 
Galaxy or nearby galaxies, are basically unex- 
plored. Such environments are particularly im- 
portant for low-mass star/BD formation studies 
because they present properties (stellar density, 
total cluster mass, number of ionizing OB stars, 
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metallicity, etc.) very different from nearby stellar 
populations and allow us to probe/discard the nu- 
merou s scenarios proposed for BD formation (e.g., 



Joergendll2005h . 



To overcome these limitations, new BD clas- 
sification methods have been devis ed on the ba- 
sis of near-IR broad-band (see, e.g., Leggett et al.l 
2OIOI : Schmidt et al. 2010l and references therein) 



and narrow-band photometry that targets unique 
molecular features of L and T dwa rfs, in particu- 
lar H2O and CH4 b ands (see, e.g., Gorlova 20031 : 



Mainzer et al.1120041 ). These techniques are equiv- 



alent to extremely low resolution spectroscopy and 
can be confidently applied for statistical purposes, 
e.g. to detect and classify bona fide very low-mass 
stars and BDs in large imaging surveys. However, 
such samples are inevitably affected by spectral 
contamination, since ground-based measurements 
of the depth of water bands are complicated by the 
presence of variable water bands in the Earth's at- 
mosphere itself. Thus, these kinds of studies are 
more efficiently conducted from space and, indeed, 
HST/NICMOS photometry has been aheady used 
to measure the strength of w a.ter bands for BP 
spectral classification purposes (jNaiita et al ]l2000t 



Andersen et al.l 120061 ). However, NIC MOS's field 



of vie w and sensitivity are still limited ([Thompson 



19981). The new Wide Field Camera 3 (WFC3) on 



board the HST is very sensitive at both optical and 
IR wavelengths and brings unique capabilities to 
star formation studies. Specifically, the WFC3 is 
equipped with a set of four near-IR narrow-band 
filters spanning the J and H bands and covering 
H2O molecular features useful for BD character- 
ization, overcoming the contamination by atmo- 
spheric water bands affecting ground-based mea- 
surements and exceeding the sensitivity and field 
coverage of NICMOS. 

In this paper we present WFC3/IR observa- 
tions, including narrow-band imaging in the H2O 
bands, of the young massive cluster in the core 
of NGC 3603. NGC 3603 (i?Aj2ooo « ll^lS'"?", 
Dec J2000 ~ — 61°15'30") is a gala ctic giant HII 
region at a distance of 6-7 kpc ( Melena et al.l 
.20081 . a nd references there in). The young (1- 
20 Mvrs. lSeccari et "al]l2010[ ) compact stellar clus- 
ter (HD97950 or NGC 3603YC) lies at the core 
of this region and has long been the center of at- 
tention for the relatively numerous population of 
massive stars (3 WNL, 6 03-type stars and nu- 
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merous late 0-type stars). Together, these re- 
gions have a bolometric luminosity of 100 times 
that of the Orion cluster and 0.1 times that of 
NGC 2070 in the 30 Doradus complex in the Large 
Magellanic Cloud and a total ma ss in excess of 



10^ M© fsee lHaravama et al.|[2008l and references 



therein). These overall characteristics make NGC 
3603 the closest small-scale resolved prototype of 
the starburst clusters commonly seen in active 
galaxies and the first natural place to probe the 
substellar population of an environment very dif- 
ferent from those observed so far. The pioneer- 
ing observations presented in this paper represent 
the "first look" at a star-burst cluster in the H2O 
bands and the first attempt to "see" into its BD 
population. 

This paper is organized as follows: Sect. [5] de- 
scribes the WFC3 observations in NGC 3603 and 
the data reduction procedure. In Sect. [3] we de- 
velop a method to identify and classify BDs on the 
basis of WFC3/IR colors in the H2O bands, dis- 
cuss its advantages and limitations in comparison 
with predictions from synthetic spectra of stars 
and BDs and use the observations in NGC 3603 
to quantify the effect of uncertainties in interstel- 
lar extinction. In Sect. 2] and [S] we apply our 
method to identify objects presenting H2O absorp- 
tion bands in the field of NGC 3603 and discuss 
the surprising properties of the selected sample. 
The conclusions of this work are given in Sect. |6l 

2. Observations and data reduction 

NGC 3603 was observed between August and 
November 2009 with the WFC3 on board the re- 
furbished HST in both the optical (UVIS) and in- 
frared (IR) channels as part of the early release 
science (ERS) program. The analysis of the UVIS 
data-set has been presented in a previous paper 
dedicated to the characterization of the stellar 
population in NGC 36 03 and to the study of its 
star formation history (jBeccari et al. 2010 ). Here 
we focus on the WFC3/IR observations. 

The WFC3 IR detector is a single 1024x1024 
HgCdTe CCD with a total field of view of 
123" X 136" at a plate scale of 0"13 per pixel. The 
observations in NGC 3603 were performed through 
the FllOW (J-band) and F160W (H-band) broad- 
band filters, the F127M, F139M and F153M 
medium-band filters, and the F128N narrow-band 



filter (see characteristics in Table [T]) . For each fil- 
ter, three images with approximately the same ex- 
posure time were taken with a constant '^3.5 pixel 
dithering in order to allow for the removal of cos- 
mic rays, hot pixels and other detector blemishes. 
The WFC3/IR detector is affected by persistence 
effects, which according to the WFC3 handbook 
are below the read noise in less than 1 5 minutes 
(see Sect. 7.9.4. bv lDressel et al.ll2010( ). In order 
to ensure that persistence spots do not compro- 
mise the photometry, especially in the center of 
the cluster where there is a high concentration of 
bright OB stars, our observations were performed 
so that the time interval between two consecutive 
visits is always longer than 1.5 day (see Table [2]). 
All observations were performed so that the core 
of NGC 3603 is roughly located at the center of 
the camera's field of view (FoV) . In Figure [1] we 
show a mosaic of the images in the FllOW filter 
as obtained with the PyRAF/MULTIDRIZZLE 
package. 

The photometric analysis of the entire data-set 
was performed on the flat-fielded (FLT) images 
multiplied by the pixel-area maps to get unifor- 
mity in the measured counts of an object across 
the field. About 150 isolated and well-exposed 
stars were selected in every image over the entire 
FoV in order to properly model the point spread 
functi on (PSF) with the DAOPHOTII/PSF rou- 
tine (|Stetsonl Il987t l. We find that a first or- 
der polynomial gives a good approximation to 
the PSF variation across the FoV. A first list of 
stars was generated by searching for objects above 
the 3(T detection limit in each individual image 
and a preliminary PSF fitting was performed us- 
ing DAOPHOTII/ALLSTAR. We then used DAO- 
MATCH and DAOMASTER to match all stars in 
each chip, regardless of the filter, in order to get an 
accurate coordinate transformation between the 
frames. A master star list was created using stars 
detected in the FllOW image (the deepest of the 
IR data-set) with the requirement that a star had 
to be detected in at least two of the three im- 
ages in this filter. We used the sharpness (sh) 
and chi square (chi) parameters given by ALL- 
STAR to remove spurious detections. The final 
catalogue was then obtained by rejecting residual 
spurious sources (~30%) through visual inspection 
of the drizzled images. The maste r list was the n 
used as an input for ALLFRAME (IStetsonlll99l . 
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Fig. 1.— WFC3 2' X 2' mosaic image of NGC 3603 in the FllOW filter as obtained with the 
PyRAF/MULTIDRIZZLE package. The open symbols indicate the positions of the 9 objects showing water- 
vapor absorption discussed in Sect. |4] and [5] (squares indicate the 3 objects whose photometry might be 
contaminated by nearby bright/saturated stars). All of them are located within a radius of ^^30 arcsec from 
the cluster center (i.e. ~1 pc at the cluster distance). Object ID numbers are as in Tabled The diffuse 
emission in the lower part of the image is due to scattered light from the reflection nebula surrounding the 
cluster and to the nebular Paschcn-/? emission line (1.2818 /im). 
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which simultaneously determines the brightness of 
the stars in all frames while enforcing one set of 
centroids and one transformation between all im- 
ages. All the magnitude values for each star were 
normalized to a reference frame and averaged to- 
gether, and the photometric error was derived as 
the standard deviation of the repeated measure- 
ments. The final catalogue contains photometry 
in the FllOW, F160W, F127M, F128N, F139M 
and F153M filters for 9693 stars. 

The photometric calibration to the standard 
VEGAM AG photometric sy stem was performed 
following Kalirai et al. ( 2003) • A sample of bright 
isolated stars was used to correct the instrumental 
magnitudes to a fixed aperture of 0'.'4. The magni- 
tudes were then transformed into the VEGAMAG 
system by adopting the synt hetic zero points for 
the WFC3/IR bands given bv lKalirai etall (|2009l) 
(see their Table 5). 

The WFC3 /IR channel is affected by geometric 
distortion and a correction is necessary in order to 
properly derive the absolute positions of individ- 
ual stars in each catalogu e. We used the distortion 
coefficients derived by iKozhurina-Platais et al 



()2009l) ;o obtain relative star positions that are 
corrected for distortion. We then used the stars 
in common between our IR and the 2MASS cat- 
alogues to derive an astrometric solution and ob- 
tain the absolute R.A. and Dec. positions of our 
stars. We find a systematic residual of ^0.3" with 
respect to the 2MASS coordinates. 

In Table [2] we summarize the exposure times, 
saturation limits and limiting magnitudes at the 
3(7 level in each filter. At the distance (6-7 kpc) 
and age (1-20 Myr) of NGC 3603 and in the 
absence of interstellar extinction, these limit- 
ing magnitudes correspond to cluster members 
with masses down to 0.15—0.2 Mq, according 
to the pre-m a in seq uence (PMS) isochrones by 
Baraffe et al.l ()l998 ) specifically calculated for 
the HST/WFC3 filterfl Considering the typi- 
cal extinction of Ay =5. 5 mag toward the center 
of NGC 3603 (Sect. [SJ]), the mass limit would 
be ~0.4 Mq. The large uncertaint y on the dis- 
tance to NGC 3603 (~1 kpc or more: lMelena et al 



20081 ) results in an uncertainty on our photometric 
limits of about 0.5 mag, i.e. '-^0.05 Mq for clus- 
ter members. This means that our observations 



probe the cluster population down to the substel- 
lar boundary in the ideal case of no extinction and 
should be too shallow to detect sub-stellar mem- 
bers. However, as we will see in Sect. 0] and [Sj we 
identified a sample of very cool objects presenting 
H2O absorption bands that are likely members of 
NGC 3603. Our method to identify objects with 
H2O in absorption is based on the use of color 
indices (Sect.[3|) and, hence, is not affected by the 
uncertainty on the cluster distance. 

3. Using WFC3 narrov^r-band photometry 
to characterize brown dwarfs 

The WFC3 is equipped with a set of four 
near-IR medium and narrow-band filters centered 
at 1.27/im (F127M), 1.28Afm (F128N), 1.39/.tm 
(F139M), and 1.53/im (F153M). The F139M fil- 
ter probes the II2O absorption band centered at 
~1.4/im while the F153M filter measures the edge 
of this band (Figure ^ and is also sensitive to 
methane absorption. Other molecular absorption 
lines (CN, HCN, metal hydrides and alkali met- 
als) characterize the spectrum of very cool objects 
(Te// <3000 K) between 1.3 and ^.5^lul, how- 
ever H2O is by far the strong est molecular ab- 



sorbe r in th is wavelength range (jKirkpatrick et al 
19991 l2005h . Because this H2O featu re is unique 



^Available at |http:/ /perso.ens-lyon.fr/isabelle.barafFe/| 



to ve ry cool objects (see Figure 1 in iMainzer et al.l 
20041 ). the F139M and F153M filters can be used 
to identify very low-mass stars and BDs and es- 
timate their approximate effective temperature 
by measuring the strength of the water absorp- 
tion bands, provided the stellar continuum in 
this wavelength region is known. This informa- 
tion is easily obtained by combining measure- 
ments in the other two filters. Indeed, F127M 
can be used as the primary continuum filter for 
it covers a wavelength range relatively featureless 
in late- type objects; the only contamination ex- 
pected in this wavelength range comes from the 
Paschen-/? line (1.2818/im), which is expected to 
be in emission in young accreting objects (see, e.g., 
Folha fc EmersonI l200l[ ) and is present in nebu - 
lous regions as well (see, e.g., iKennicutO Il998l ). 
The narrow F128N filter, centered at 1.28/xm with 
a bandwidth of 15.9 nm, is a good indicator of 
the possible contamination due to the Paschen- 
13 emission line. Thus, the flux ratio between the 
F139M and F153W bands with respect to the con- 
tinuum, corrected for Paschen-/3 emission, gives a 
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F187M F128N F139M 



Fig. 2.— The F127M, F128N, F139M and 
F153M WFC3/IR transmission bands. An ex- 
ample of a normalized AMES-Co nd model BD 
spectrum ( Hauschildt et al. 19991 ) with Toff = 
2000 K and logg = 4.0 is over-plotted (black 
line). The broad H2O absorption band centered 
at ~1.4/xm is clearly visible. For comparison, 
we plot the AMES-Con d model stellar spectrum 
(|Hauschildt et al.lll999 ') with Toff = 3500 K (grey 
line), which does not show the H2O feature. 



Table 1: Characteristics of the WFC3/IR fihers 
used in this paper: pivot wavelength (Ap), band- 
width (AA) and photometric zero-point for the cal- 
ibration to the VEGAMAG standard system. 



Filter 


Ap^ 


AA^ 


Zp*^ 




(nm) 


(nm) 




FllOW 


1153.4 


443.0 


26.07 


F127M 


1274.0 


68.8 


23.69 


F128N 


1283.2 


15.9 


21.98 


F139M 


1383.8 


64.3 


23.42 


F153M 


1532.2 


68.5 


23.21 


F160W 


1536.9 


268.3 


24.70 



"As defined in Sect. 7.5 of lDressel et al.l 1I2OIOI '). 
''From .Kalirai et al.l (120091) . 



good measure of the strength of the water absorp- 
tion bands. 

Our images are produced in units of elec- 
trons/sec and, hence, the counts scale with the 
bandwidth of the filter. Thus, to compute the 
flux density (Flux) in the continuum and the wa- 
ter line, we divided the corresponding counts (C) 
by the effective bandwidths of the filters: 



cont = 



Cpi 



27M 



28N 



AApi27M ~ AApi287V 
Cf139M 



FluXF139M = 



FIuXfi53M = 



AApi39M 

Cf153M 
^^F153M 



(1) 



(2) 



(3) 



where AApi27M, AApi287v, AXfi39M and 
AApi53j\/ are the bandwidths of the four filters 
reported in Table [1] Equations 1-3 implicitly as- 
sume that the adopted filters have a rectangle 
response curve (top- hat filters), which is realistic 
in the case of WFC3/IR filters (Figure [l]). 

We then convert the flux density in magni- 
tude (mag) and define two narrow-band indices, 
l{lS9)=magcont-magFi39M and I(153)=magcont- 
TnagFi53M, which are a measure of the water ab- 
sorption and can then be used to identify very 
late-type objects and estimate their effective tem- 
perature. 

In order to verify this approach, we built an 
empirical relation between these color indices and 
the effective temperature by using the actual spec- 
tra of 552 stars and BDs available from the liter- 
ature. These objects range in spectral type from 
05 to T9 (Te// =45000-600 K) and were drawn 
from the Bruzual-Persson-Gunn-Stryker (BPGS ) 
atlas (jStrecker et al.lll979HGunn fc Strvkeij|l983[) . 



the Brown Dwarfs Spectroscopic Survey (BDSS, 
McLean et al1|2003l ), the SpeX Prism Spectral Li- 
braries (SpeXjl and the spectroscopi c databases 
for ve r y low-mass s tars and BDs by iReid et all 
(|200l[ ). iTestil (|2009i) . and B. Burningham (pri- 
vate communication). Using the WFC3/IR fil- 
ter transmission curves and the STSDAS/Synphot 
synthetic photometry package in IRAF, we calcu- 
lated the magnitudes of each object in each band 



^http;/ /www. browndwarfs.org/spexprism/ 
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by integrating the total flux in that band. Mag- 
nitudes are given in the VEGAMAG photomet- 
ric system using the zero points for the WFC3/IR 
filters listed in Table [TJ Figure [3] shows the re- 
lation between the two water-band indices, de- 
rived from the observed spectra, and the effec- 
tive temperatures (Te//) of the 137 stars and 
BDs in our sample. The figure shows that there 
is reasonable agreement between all of the sets 
of observations and the thick solid line repre- 
sents the mean relation derived by averaging the 
measurements within steps of 250 K. Down to 
Te// ~3500 K, both the 1(139) and 1(153) color 
indices have a mean value close to zero, indicating 
that no water is present in the spectrum. The on- 
set of the water- vapor absorption feature is clearly 
seen below Te// ~3500K, where both color in- 
dices become more and more negative as Te/ / de- 
creases. This decrease is roughly monotonic in 
the range 3500— 500K for 1(139) and in the range 
2500— 500K for 1(153) and, hence, these color in- 
dices can be used for BD spectral classification 
purposes. 

We compared this empirical relation with 
expectations from model spectra of very-low 
mass stars and BDs. This can give further in- 
sight into the physical processes occurring in 
their atmosphere, w hich are still not clear (see 
Ghabrier et al-llioOOl and references therein). We 



used the NextGen, AMES-Dusty and AMES-Gond 
model spectra which are, to our knowledge, the 
most complete grids of model atmospheres for late- 
type objects currently available. The NextGen 
models are dust-free model atmospheres (see 
iHauschi ldt et all Il999l). The AM ES-Dustv and 
AMES-Cond models (jAllard et al.l lioOl) take into 
account the formation of condensed species, which 
depletes the gas phase of a number of molecular 
species and refractory elements, modifying signif- 
icantly the atmospheric structure and, hence, the 
emergent spectrum. While in the AMES-Dusty 
models the condensed species are included both 
in the equation of state and in the opacity, tak- 
ing into account dust scattering and absorption, 
in the AMES-Cond models the opacity of these 
condensates is ignored, in order to mimic a rapid 
gravitational settling of all grains below the pho- 
tosphere. For furth er details about t hese grids of 
models we refer to Ghabrier et al. (|2000l ). The 
color indices of the model spectra were derived 



by folding them with the WFG3/IR filters, as was 
done with the actual spectra. The grey-filled areas 
in Figure |3] show the color vs. Te// dependency 
for the three sets of models. The comparison be- 
tween the model predictions and the spectroscopic 
data shows that: 

1. The relation between both our narrow-band 
color indices and Te// depends on stellar 
gravity, as indicated by the shaded areas in 
Figure [3l which encompasses all the calibra- 
tions for models in the range 3.5 < log 5 < 
5.5. This variation with logg is more signif- 
icant for the 1(153) color index, but in any 
case is well within the spread of the spectro- 
scopic data; 

2. the predictions by the three sets of models 
are in good agreement with the actual spec- 
tra down to -3GGGK; 

3. below ~3000K the predictions of the three 
sets of models differ significantly from one 
another because the effect of dust grain for- 
mation in the atmospheres is treated differ- 
ently, as explained above, and below 2000K 
they clearly disagree with the actual spec- 
tra. Data from the actual spectra show a 
linear decrease of both color indices with 
Te//, more similar to the prediction by the 
AMES-Gond models. However, for both 
color indices this decrease is less steep than 
that predicted by the AMES-Gond models, 
mimicking an average behavior between the 
AMES-Gond and the AMES-Dusty models. 
In other words, the data support a scenario 
where both dust gravitational settling and 
scattering/absorption operate in the atmo- 
sphere of very cool objects and, hence, affect 
their spectrum. 



AUard et al.l (|200ll) suggest the use of the 
NextGen models for Te// >2700K, the AMES- 
Dusty models for 1700K< Te// <2700K, and the 
AMES-Gond models for Te// <1400K. Between 
1700 and 1400K, a new set of models (AMES- 
Settl) are being studied and are still in the ex- 
perimental phas^. Indeed, Figure [3] shows that 



3 See R AUard 

|http://perso.ens- lyon.fr/fr ance.allard/l 



web-page: 
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Fig. 3. — 1(139) and 1(153) color indices as a function of the effective temperature for stars and BDs 
with actua l spec tra from the BPGS atla s (stars), the BDSS (circles), the SpeX libraries (plus symbols), 
Reid et all ( 200lh (triangles) . iTesti ( 2009t ) (upside down triangles) and B. Burningham (private communica- 
tion, squares). The solid line is the as-measured calibration relation reported in Table[3l The grey- filled areas 
represent the relations derived from the m odel spectra, as indicat ed in the legend. The thin and thick dashed 
lines represent the calibration relation bv lAaronson et ah (1978) for giant and dwarf stars, respectively. 
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this prescription roughly reproduces the trend ob- 
tained from the observed spectra. In particular, 
the AMES-Dusty models are roughly consistent 
with the 1(139) and 1(153) colors for 1700K< 
Teff <2000K, while the AMES-Cond models un- 
derestimate both colors for Teff <1400K for any 
stellar gravity in the range 3.5 < \ogg < 5.5. 
Thus, given the large uncertainties affecting model 
predictions in this temperature range, we conclude 
that the mean relation between the water-band in- 
dices and Te/ / derived from the observed spectra 
(solid line in Figure [3|) is more reliable for the 
identification and classification of BDs. We re- 
port this calibration relation in tabular form in 
Table [31 the uncertainties on the color indices are 
calculated as the dispersion of the data in each 
temperature bin. 

We also show in Figu re El the ground-based 
H2O indices computed by Aaronson et al. ( 19781 ) 
as a function of effective temperature for giant and 
dwarf stars. These indices were used for quite 
some time in constructing population synthesis 
models for galaxies and clusters. It is then re- 
assuring to see that our 1(139) color index per- 
fectly reproduces the trend of the lAaronson et al.l 
(|l978 l) II2O indices, while the 1(153) only slightly 
underestimates them. Note that this agreement 
is not straightforward because the measurements 
were taken from different facilities (one ground- 
based, the other a satellite), different filters, dif- 
ferent continuum estimates, etc. 

Figure |4] shows the 1(139) vs. 1(153) color-color 
plot derived from our calibration (thick solid line) . 
This diagram shows that, when no information on 
the effective temperature of an object is available, 
WFC3 photometry in the four near-IR bands can 
be used to give an estimate of its Teff value. In 
particular, objects close to or below the sub-stellar 
limit {Teff ~3000K) are expected to be located 
in a region of this diagram well separated from 
that of more massive stars. We also used our BD 
reference sample to estimate the accuracy on the 
color-derived Te//. We performed a linear fit of 
the spectroscopic Te// of the BDs in our refer- 
ence sample as a function of the Te// predicted 
by the color— color relation. The computed RMS 
indicates that the two sets of measurements agree 
within ^^350 K. Thus, the accuracy of our method 
is sufRcient for BD identification purposes, even 
though it provides only a rough estimate of their 



Table 2: Start time, total exposure time, satu- 
ration limit and limiting magnitude at 3(7 level 
in each WFC3/IR filter for the observations in 
NGC 3603. 



Filter 


Dithering 


Star time 
(Julian day) 


Tex 

(sec) 


Mag Sat. 


Mag 3a 


FllOW 


1 
2 
3 


55070.226 
55070.291 
55070.357 


3x200 


12.0 


20.5 


F127M 


1 

2 
3 


55070.241 
55070.306 
55070.372 


3x800 


10.0 


20.0 


F128N 


1 
2 
3 


55070.251 
55070.315 
55070.382 


3x400 


9.0 


19.5 


F139M 


1 
2 
3 


55070.453 
55070.519 
55070.586 


3x800 


10.0 


20.0 


F153M 


1 
2 
3 


55070.424 
55070.491 
55070.557 


3x800 


9.5 


19.5 


F160W 


1 
2 
3 


55070.229 
55070.294 
55070.360 


3x200 


11.5 


19.5 



Table 3: As-measured calibration relation between 
the water-band indices 1(139) and 1(153) and the 
effective temperature. 



Spec.^ 




1(139) 


1(153) 


Type 


(K) 






T9 


500 


-5.70±0.34 


-1.68±0.05 


T8 


750 


-4.43±0.38 


-1.27±0.01 


T6 


1000 


-3.17±0.41 


-0.85±0.08 


T4.5 


1250 


-1.90±0.45 


-0.43±0.15 


L6.5 


1500 


-1.08±0.33 


-0.19±0.14 


L4.5 


1750 


-0.56±0.11 


-0.05±0.08 


L2.5 


2000 


-0.49±0.09 


-0.08±0.07 


LI. 5 


2250 


-0.44±0.10 


-0.13±0.07 


LO 


2500 


-0.39±0.10 


-0.18±0.06 


M8 


2750 


-0.33±0.10 


-0.23±0.06 


M6 


3000 


-0.21±0.09 


-0.17±0.06 


M4 


3250 


-0.10±0.08 


-0.10±0.06 


M2.5 


3500 


-0.04±0.05 


-0.05±0.03 


MO. 5 


3750 


0.01±0.02 


-O.OliO.Ol 


K7.5 


4000 


O.OlifcO.Ol 


-0.02±0.01 



^ Spectral type s are inferr e d usin g the temper 
ature scal e bv IVrba et al. 



dwarfs, Luhman et al.l 



Kenyon fc Hartmann |l995f ) for K dwarfs. 



( 20041 ) for L and T 
j200j) for M dwarfs and 



9 



effective temperature. Two caveats must be taken 
into account: i) our Teff estimate is based on the 
assumption that uncertainties on magnitude are 
not higher than 10% in each of the four filters; ii) 
there has been a controversy in recent years over 
the spectroscopic effective temperatures scale for 
L and T dwarfs in that differences of the order of 
250 K ex ist amo ng different cal ibrations (see, i.e., 



Basri et a l. 2000: ; Schweitzer et al. 2001. 120021 ) 



Since our calibration relation is based on the use 
of near-lR filters that are very close to each other 
in central wavelength, it is only slightly affected by 
uncertainties on interstellar extinction (Figure HJ. 
This makes our method particularly suitable to 
identifying BD candidates in young clusters, where 
extinction is normally highly variable. This issue 
is further discussed in Sect. 



3.1. 



Using the J-band as the primary con- 
tinuum filter 



In this section we discuss the adoption of the 
WFC3 /IR J-like bands as the primary continuum 
indicators when computing the 1(139) and 1(153) 
color indices. 

As explained in Sect.[3l the use of the F128N fil- 
ter allows us to correct the continuum emission for 
the possible contamination due to the Paschen-/? 
line, which is expected to be in emission in young 



accre ting objects (see, e.g., iFolha fc Emerson 
20011) . However, since F128N is a narrow-band 
filter (AA =13.54 nm), long exposure times are 
required to reach a photometric depth suitable 
for BD detection in distant star forming regions. 
Using the WFC3 /IR-imaging exposure time calcu- 
latoilfl, we estimate that observations through the 
narrow F128N filter typically require an exposure 
time five times longer than for the three medium- 
band filters (F127M, F139M and F153M) to reach 
the same photometric depth. Thus, observations 
through the F128N filter are very time consuming 
and dominate the required observing time. 

For those objects which are not expected to ac- 
crete (e.g. field BDs) and are not located in nebu- 
lar regions, the J-band filter can be used as a con- 
tinuum filter. The WFC3/IR camera is equipped 
with two J-like fihers: FllOW and F125W. In Ta- 
ble |4] we provide the calibration relations between 



*http: / /etc. stsci.edu/webetc/mainPages/wfc3IRImagingETC.jsp 



the water-band indices 1(139) and 1(153) and the 
effective temperature when the continuum is es- 
timated using these two J-like filters. Note that, 
because both FllOW and F125W are wide-band 
filters, the continuum estimates based on the use 
of these filters are slightly more sensitive to uncer- 
tainties in interstellar reddening. 

3.2. Quantifying the efftect of interstellar 
extinction 

In this section we use the observations in 
NGC 3603 to quantify the effect of uncertainties 
in interstellar extinction on our BD classification 
method. 

Figure[5]shows the dereddened 1(139) vs. 1(153) 
color-color diagram for point-like sources in the 
field of NCG 3603. In order to correct our mag- 
nitudes for interstellar extinction, we followed the 
study of differential reddening in NGC 3603 by 
Sung fc Bessell ()2004f ). These authors were able 
to map the variation of E(B-V) as a function of 
the distance from the cluster centre (see their Fig- 
ure 5b) using multi-band HST photometry of the 
bright massive stars (which are saturated in our 
images). They found that the value Ay =4.5 is rep- 
resentative of the very centre of the OB stars asso- 
ciation, while they noticed an increase toward the 
external regions. Adopting the total-to-selective 
extinction rati o (Ry) value of 3.55, suggested by 
Sung fc Bessell (|2004l ). we estimate the value of 
Ay =5.5 (i.e. Aj «1.5) to be representative of 
the mean visual extinction in the area sampled 
by our observations (from ~10" to ~70"). We 
ad opted this value of A y and the extinction law 
bv lCardein et alj I (|l989l) to correct the magnitude 
of stars in our catalogue. As shown in Figure [5l 
the dereddened photometry matches rather well 
our calibration relation, thus confirming the accu- 
racy of our method. 

As mentioned in Sect. |31 since our calibration 
relation is based on the use of near-IR filters that 
are very close to each other in central wavelength, 
it is only slightly affected b y uncertainties on in - 
terstellar reddening (see also lMainzer et aL I l2004h . 
The arrow in Figure [S] represents the reddening 
vector for Ay=l, which is the typical uncertainty 
in NGC 3603. This uncertainty in interstellar ex- 
tinction translates the data points by 0.03 mag 
in 1(139) and 0.07 mag in 1(153), i.e. within the 
photometric uncertainties. 
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1(139) 

Fig. 4. — 1(139) vs. 1(153) color-color diagram for stars and BDs with actual spectra (symbols are as in 
Figure |3]). The solid line is the as-measured calibration relation between the two color indexes reported in 
Table [31 The arrow represents the =1 reddening vector. 



11 



Table 4: As-measured calibration relation between the water-band indices 1(139) and 1(153) and the effective 
temperature when using the FllOW and F125W filters as continuum estimators. 



bpec. ' 
Type 




Continuum filter FllOW 
Hioy ) 


Hloo ) 


Continuum filter F125W 
i(loy ) 


Hioo ) 


19 


oUU 


-5.ZD±0.5o 


-i.ZO±U.ZO 


-b.Z6±[).b ( 


1 OO _l_ A O /I 


io 


YOU 


-4.1U±0.54 


A no_i_A OO 


A nO_l_A CO 

-4.(Jo±U.oo 


A AO_l_A OO 

-(j.yz±(j.zz 


io 


1000 


-2.94±0.51 


A /^O 1 A 1 


o A /I 1 A rr A 

-z.94±0.50 


A £?0 1 A OA 

-U.oz±(J.zU 






1 77-l-n 47 


A Ql ^0 90 
-u.oizcu.zu 


1 7Q+n A(\ 
- 1 . ( y inu.^u 


A QO-Ufl 1 Q 


L6.5 


1500 


-1.01±0.33 


-0.11±0.16 


-1.04±0.32 


-0.14±0.14 


L4.5 


1750 


-0.50±0.11 


0.02±0.10 


-0.53±0.10 


-0.02±0.08 


L2.5 


2000 


-0.39±0.11 


-0.01±0.11 


-0.44±0.09 


-0.05±0.08 


L1.5 


2250 


-0.34±0.14 


-0.03±0.13 


-0.41±0.11 


-0.10±0.08 


LO 


2500 


-0.32±0.13 


-O.lliO.lO 


-0.37±0.10 


-0.16±0.07 


M8 


2750 


-0.30±0.11 


-0.20±0.08 


-0.32±0.10 


-0.22±0.06 


M6 


3000 


-0.21±0.09 


-0.17±0.07 


-0.22±0.09 


-0.18±0.06 


M4 


3250 


-0.11±0.08 


-0.11±0.07 


-0.12±0.07 


-0.12±0.06 


M2.5 


3500 


-0.09±0.05 


-0.10±0.05 


-0.09±0.04 


-0.10±0.04 


M0.5 


3750 


-0.07±0.03 


-0.10±0.03 


-0.06±0.01 


-0.08±0.02 


K7.5 


4000 


-0.14±0.02 


-0.17±0.02 


-0.09±0.01 


-O.lliO.Ol 



^ Spectral types are inferred using t he temperature scale bv IVrba et al.l l|zui 
Luhman et alj ( 20031) for M dwarfs and Kenvon fc Hartmann ( 19951) for K dwarfs 



(|2004[ ) for L and T dwarfs, 



Note also that, for large values of Ay, de- 
pends on the intrinsic color of the star because of 
flux weigh ting through t he different filter response 
functions ( Wildev 1963L and references therein). 
However, the fact that we compute colors from fil- 
ters very close to each other in central wavelength 
makes this dependency moot. 

4. Selection of objects with water-vapor 
absorption in NGC 3603: what are 
they? 

As seen in Sect.jU our observations in NGC 3603 
probe the cluster population down to the substel- 
lar boundary in the best case, i.e. in the absence of 
extinction, and hence we do not expect to detect 
any object presenting H2O absorption bands be- 
longing to this cluster. However, the dereddened 
1(139) vs. 1(153) color-color diagram of the point- 
like sources in the field of NCG 3603 (Figure [5]) 
surprisingly shows a number of objects located in 
the region expected for BDs with significant ab- 
sorption in the water bands (i.e. T^efj ^3000K). 

In order to systematically isolate this group of 
points, we first considered the calibration relation 



1(139) vs. 1(153) in the substellar domain and cal- 
culated the dispersion (a spec) of the spectroscopic 
data used to derive this relation. We then selected 
all objects below T^eff =3000K having colors con- 
sistent with the calibration relation within laspec 
and whose photometric uncertainties still place 
the object within Zdspec of the mean calibration 
relation. These criteria take into account that, 
when measuring spectral line intensities, narrow- 
band photometry is necessarily less accurate than 
spectroscopy. We identified in this way 9 ob- 
jects with absorption in the water bands, shown as 
filled circles / squares in Figure [5] In Figure [6] we 
plot the position of these 9 objects on the dered- 
dened magcont vs. magFi27M ~ magpisgM color- 
magnitude diagram. There are 9 objects which are 
singled out as interesting sources because they ex- 
hibit more negative rnagFi27M — magpi^gM colors 
with respect to the remaining sources in the ob- 
served field, due to the 1.4/im absorption band. 
Note that there are a number of objects with 
'magpi27M — "magFisgM colors similar to these 9 
objects which were not selected on the basis of 
Figure [5j because they have larger photometric 
uncertainties on the 1(139) and 1(153) indices and 
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T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




I I I I I I u I I I I I i I I I I I 

-1.5 -1.0 -0.5 0.0 0.5 

1(139) 

Fig. 5. — Dereddened 1(139) vs. 1(153) color-color diagram for point-like sources in the field of NCG 3603 
(small dots) . The solid line is the calibration relation derived by us using stars and BDs with actual spectra; 
these spectroscopic data are represented as open symbols (as in Figure and the dashed lines represent their 
la dispersion in the substellar domain {Teff ^3000K). The larger filled symbols with error bars indicate the 
positions of the 9 objects showing water- vapor absorption (squares indicate the 3 objects whose photometry 
might be contaminated by nearby bright/saturated stars). The arrow represents the Ay=l reddening vector. 
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Fig. 6. — Dereddened magcont vs. magpi27M — nT-agpugM color-magnitude diagram for point-like sources 
in the field of NCG 3603 (small dots). Average photometric uncertainties are reported on the right in 
steps of 0.5 mag. The bigger symbols with error bars indicate the positions of the 9 objects showing 
water-vapor absorption (squares indicate the 3 objects whose photometry might be contaminated by nearby 
bright/saturated stars). 
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do not satisfy our selection criteria. 

In Table [5] we list the observed WFC3/IR 
photometry for the 9 objects with water- vapor 
absorption and their derived effective tempera- 
tures, which lie in the range 1700-2200K. Their 
J-band observed magnitudes lie in the range 19.5- 
21.5; assuming the distance modulus of NGC 3603 
(~14 mag) and the average interstellar reddening 
in its direction (Aj «1.5), this range corresponds 
to J-band absolute magnitudes between 4.1 and 
6.1 and to stellar luminosities in the range 0.6- 
4L0. 

Since the photometric sensitivity of our data 
does not allow us to detect NGC 3603 members 
below 0.15-0.2 M0 in the ideal case of no extinc- 
tion, we now investigate what these objects are. 
Before exploring the various possibilities, we col- 
lect further hints on the nature of these objects 
making use of our observations in the broad-band 
FllOW (J-band) and F160W (H-band) filters: 

1. Figure[7]shows the dereddened J vs. {J~H) 
color-magnitude diagram (CMD) for the 
point-like objects detected in NGC 3603 
together wit h the zero age main sequence 
(ZAMS) of iMarigo et all (|2008[ ). which 
roughly separates the locus expected for 
field main-sequence (MS) stars from the 
PMS star locus, and th e PMS isochrones and 
evol utionary tr a cks b y iBarafFe et al. I (Il998h 
and Siess et al. (|2000[ ). which take H2O and 
other prominent near-IR absorption features 
into accouniH. The position of the 9 objects 
with water absorption in this CMD indi- 
cates that their luminosity is comparable 
to that of PMS/MS stars with masses be- 
tween ^0.2 — 1 Mq, i.e. much higher than 
expected for PMS objects with a surface 
temperature in the range 1700-2200K. More- 
over, only a few of them have (J — H) col- 
ors typical of such cool PMS objects, while 
the majority appear to be bluer and, hence, 
older. While the over-luminosity of these 9 
objects has no straightforward explanation 
(see Sect. [S]), their pecuhar {J — H) colors 
can be partly explained considering the star 



®Note that isochrones have been corrected for a constant 
color term of J — _ff=0.18 mag to account for the difference 
between the standard JH Johnson's filter and the WFC3 
FllOW and F160W filters. 



formation h istory of NGC 3603. Indeed, the 
analysis by Beccari et al. ( 2010l ) focused on 
stars with Ha excess emission, a robust in - 
dicator of youth (e.g.. White &: Basri 2003[) . 
found that star formation in and around the 
cluster has been ongoing for at least 10-20 
Myr; 2/3 of the young population have ages 
from 1 to 10 Myr with a median value of 3 
Myr according to their position on the CMD, 
while 1/3 of them lie at or near the ZAMS, 
thus suggesting a considerably higher age, 
of order of 20-30 Myr. The 9 objects with 
water-vapor absorption appear to follow the 
same behavior, being spread between the 
PMS region of the CMD and the ZAMS. Al- 
though the statistic of this sample is hardly 
significant, we note that percentages are 
reversed with respect to the Ha emitters 
population, i.e. '^2/3 of the objects with 
absorption in the water bands lie near the 
ZAMS and only ~l/3 (i.e. 2-3 objects) are 
located in the PMS region; 

The spatial distribution of the 9 objects 
with water-vapor absorption appears to be 
strongly clustered towards the luminous core 
of NGC 3603 (see Figure [l]). We attempt 
to use a Kolmogorov-Smirnov (KS) test 
to check whether the radial distribution of 
these 9 objects with respect to the cluster 
center is consi stent with that of t he cluster 
members (see Beccari et al.l[201ol for a de- 
tailed description of this test). The small 
number statistic makes the results of this 
test not highly reliable; however, the general 
indication is that these 9 objects are more 
centrally concentrated than the field pop- 
ulation, mimicking the radial distribution 
of the young stars belonging to the clus- 
ter. We also noticed that all 9 objects are 
preferentially located on the North side of 
NGC 3603 (Figure [U. The presence of pil- 
lars in the South region suggests a higher 
gas and dust density which could prevent 
us from detecting faint objects on this side 
of the cluster. A high resolution extinc- 
tion map of NGC 3603 has been recently 
computed (X. Pang & A. Pasquali, private 
communication) us ing WFC3 narrow- band 
imaging in the Ha ( Beccari et al. 2010l) and 
Pa/3 lines (this work). The very preliminary 
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Table 5: Observed WFC3/IR photometry and estimated elFective temperature for the 9 objects showing 
water-vapor absorption. 



ID 


RA .J2000) 


DEC J2000 




"^F127A/ 


"^F128Af 


"^F139A-/ 


"^F153A-/ 


"^F160W 




Spectral 




(hh:mm:=s) 


(dd:mm:sB) 














(K) 


Type 


1 


11:15:04.79 


-61:15:42.85 


20.65±0.12 


20.27±0.11 


20.1S±0.15 


20.43±0.10 


19.66±0.11 


19.58±0.09 


2100 


L2 


2 


11:15:04.93 


-61:15:42.75 


20.43±0.13 


19.98±0.14 


20.08±0.17 


20.13±0.0S 


19.35±0.11 


19.23±0.09 


2000 


L2.5 


3 


11:15:06.17 


-61:15:17.76 


21.44±0.06 


20.83±0.08 


21.13±0.40 


20.86±0.11 


20.20±0.08 


20.22±0.07 


2200 


LI. 5 


4* 


11:15:00.82 


-61:15:19.26 


19.73±0.23 


19.48±0.14 


19.55±0.14 


19.86±0.16 


18.85±0.09 


18.67±0.12 


1700 


L5 


5 


11:15:06.98 


-61:15:26.49 


19.91±0.11 


19.37±0.09 


19.54±0.12 


19.74±0.12 


18.45±0.10 


18.32±0.10 


1700 


L5 


6 


11:15:07.00 


-61:15:25.98 


20.46±0.15 


19.66±0.14 


19.64±0.17 


19.8S±0.14 


19.05±0.09 


18.79±0.09 


1900 


L3.5 


7 


11:15:07.06 


-61:14:59.93 


20.74±0.07 


20.05±0.09 


20.02±0.13 


20.19±0.06 


19.35±0.0S 


19.36±0.06 


2100 


L2 


8* 


11:15:08.04 


-61:15:29.59 


19.63±0.11 


19.36±0.10 


19.51±0.17 


19.47±0.16 


18.61±0.10 


18.41±0.09 


2100 


L2 


9* 


11:15:09.85 


-61:15:34.44 


21.46±0.15 


20.70±0.1S 


20.49±0.25 


21.08±0.28 


19.89±0.11 


19.92±0.11 


1700 


L5 



Spectral types are inferred using the temperature scale bv IVrba et alj ()20Q4i r 



Possible photometric contamination from nearby bright/saturated stars. 



result of this extinction study does not show 
a significant difference in reddening between 
the north and south regions of the cluster 
for the depth probed by the Ha and Pa/3 
lines. However, the higher gas density in the 
South region can shield circumstellar disks 
from the radiation field of NGC 3603 more 
efhciently. This implies that circumstellar 
disks in the North and South region of the 
cluster might undergo different evolutionary 
paths. Indeed, as we well see in Sect. [SJ the 
observed properties of our 9 objects may be 
related to a specific disk dispersal mecha- 
nism, i.e. planet formation. 

Further clues as to the evolutionary phase of 
these 9 objects were obtained from optical broad- 
band and narrow-band Ha photometry, which 
would allow us to put constraints on the mass 
accretion p rocess typical of very young objects 
(see, e.g., iDe Marchi et al. 2010l ). We cross- 
matched our WFC3/IR catalog with the optical 
WFC3/UVIS catalo g of sources in NCG 3603 by 
Beccari et all (|2010l ). However, the 9 objects in 
exam are either not detected in the optical or do 
not have suffici ent photometric accura, cy to apply 
the method by De Marchi et al. ( 2010[ ). which re- 
quires detection in the three V, I and Ha bands to 
determine the mass accretion rate of PMS objects. 

We also searched for previous ground- and 
space-based observations of these 9 objects that 
could provide some hint on their variability. 

Two ground-based imaging surveys of NCG 3603 
are available to date: the J HK survey conducted 
with IS AAC@ESO-VLT by iBrandl et all (|l999h 
(see also A ndersen 20051 ) and the deep, high an- 
gular resolution JHKL' survey conducted with 



NA0S -C0N1CA@ES0-VLT by iHaravama eFal 
(j2008l ). The spatial resolution and sensitivity of 
the ISAAC data is 2-3 times poorer than for the 
WFC3 data. Thus, we could not retrieve accurate 
photometry of our 9 objects, which are located 
very close to the cluster center. However, all 9 
objects are visible, though very faint, in K-hand 
images. In the J and H-hand images, only 3 
of them are detected with sufficient accuracy to 
attempt photometry and the magnitudes we re- 
trieved are consistent with the WFC3 ones within 
the photometric uncertainties. Only 3 of our 9 
objects (namely ID 5, ID 6 and ID 8) are lo- 
cated in the central 28"x28" area of NGC 3603 
observed with NAOS-CONICA. Moreover, these 
observations have a detection limit of Jw20 (see 
Figure 3-4 in lHaravama et al.l 120081 ) . while our 9 
objects have magnitudes very close to this limit 
or fainter. Thus, the ground-based observations 
available to date for NGC 3603 do not complement 
our WFC3-IR dataset. 

NGC 3603 has been observed both with ACS 
(Advanced Camera for Surveys) and NICMOS 
(Near Infrared Camera and Multi-Object Spec- 
trometer) on board of the HST using several nar- 
row and wide-band filters spanning the 0.22-1.6/zm 
wavelength range. These observations were per- 
formed on 2005 December 29 and 1997 Septem- 
ber 24, respectively. We retrieved the multi- 
drizzled images of this dataset from the HST 
Legacy Archive^. Our 9 objects are not detected 
in the NICMOS FllOW and F160W images be- 
cause of their limit ed photometric depth (see also 



Moffat et al.l 120041 ) ■ but are detected in the ACS 



images in the F435W and F850LP filters and their 



^http: / /hla. stsci.edu/ 
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magnitudes are reported in Table |6l 

(|20inf ) 



Recently, iRochau et al 



presented a 

proper motion study of the core of NGC 3603 
based on HST/Wide Field Planetary Camera 2 
(WFPC2) observations obtained 10 years apart, 
respectively, in 1997 and 2007. Only 3 of the 9 
objects with water absorption fall in the 34"x34" 
area observed with the Planetary Camera (PC) 
and none of them is detected, most likely because 
their visual magnitude is fa inter than the detec- 
tion limit (m555 <25) of the Rochau et al. ( 2010l ) 
dataset. 

5. Discussion 

Our exploratory narrow-band survey of NGC 3603 
has led to the surprising identification of 9 ob- 
jects considerably more luminous than BDs which 
nonetheless display water-vapor absorption bands 
typical of BDs. Here we discuss possible explana- 
tions: 



1. Contamination from field objects: The 

most immediate explanation would be that 
the 9 objects presenting water- vapor absorp- 
tion are field objects unrelated to NGC 3603. 
Under this hypothesis, these 9 objects 
could be either foreground field BDs or 
background early M-type super-giant stars, 
where the pr esence of H2 O absorption was 
confirmed by iTsuiil (|2000t) . Moreover, long 
period variables and late M-type super- 
giants in metal poor environments can have 
strong H2O absorption and, hence, bluer 
J — H col ors than expected for their T^ff 
(see, e.g., Elias et al] Il985l : iLebzelter et al. 
2010t ). All these hypotheses would explain 



Table 6: Observed ACS/Wide Field Camera pho- 
tometry for the 9 objects showing water- vapor ab- 
sorption. 



ID 



niij'saoz.p 



25.863±0.797 

27.597±1.839 
26.659±1.171 
28.229±2.453 
25.532±0.728 



23.603±0.547 
23.243±0.465 
23.625±0.555 
22.437±0.323 
21.685±0.229 
21.820±0.244 
23.426±0.505 
22.036±0.282 
23.772±0.601 



why these objects appear over-luminous 
with respect to cluster BDs and why they 
are mainly located on the ZAMS in the 
CMD (Figure [T]), but on the other hand 
there are a number of objections to it: a) 
there is no explanation for the concentra- 
tion of the contaminants around the center 
of NGC 3603 (Figured]); b) metallicity stud- 
ies in the direction of NGC 3603 confirm 
that the typical Galactic solar metallicity 
(z = 0.020) i s appr o priate for t his region 



( Smartt et aD l2002t iPeimbert et al. 200 
Lebouteiller et al T l2008[ ): c^ th ^bin et 
(l2003h model of the Galactic stellar distri 



io3 
ID 



bution does not predict any M-type super- 
giants in the 0.002 deg^ area observed to- 
wards NGC 3603; d) usi ng the magnitude 
limits given in TableEl the lCruz et alT(|2007l ) 
luminosity function, and the dis tribution of 
L and T dwarfs in the Galaxy bv lRvan et al 



()2005l) . we find a predicted count in our field 
of 0.16 L-dwarfs for no extinction and 0.04 
if the line-of-sight extinction is Ay=5.5 mag 
(Sect. l3^ . However, Galactic counts mod- 
els become more and more inaccurate to- 
wards the very low-mass regime and the 
predicted number depends on the adopted 
model. On the other hand, using other 
models (|D'Antonal[l999l : iRobin et al.ll2003h 
would also result in no more than 1 or 2 LT 
dwarfs in this field. Thus, we do not expect 
to pick up any appreciable number of inter- 
loping field LT-dwarfs in the observed area 
towards NGC 3603. For all these reasons. 
Galactic field stars do not appear to be a 
major source of contamination of our sam- 
ple of objects with water-vapor absorption. 

Extra-galactic contamination: It is gen- 
erally very hard to distinguish between an 
extended sources and background compact 
galaxies on the ba sis of broad-ban d col- 



ors alone (see, e.g., Yan et al. 2008[ ). Our 



sample of water- vapor absorbers might then 
be contaminated by extra-galactic objects. 
From the most recent J-band galaxy num- 
ber counts__(see_JI^Me_J__^aiid_Figure 17 
by ICristobal-Hornillos et al. 20091 ) , the ex- 
pected number of background galaxies with 
19. 5< J <21.5, i.e. the magnitude range of 
the 9 objects under examination, in the ob- 



17 




Fig. 7. — Dereddened J vs. (J — H) color-magnitude diagram for point-like sources in the field of NCG 3603 
(small dots). Avera ge photometry uncer tainties are reported on the left in steps of 0.5 mag . The thicker solid 



hue is the ZAMS oflMarigo et al 



y uncer 



for masses M>1.4 Mq and lBaraffe et al 
the evolutionary tracks by 



while the dashed lines are the PMS isochrones of lSiess et al, 
998h for masses M<1.4 



(l200d) 



see legend) . The thiner solid lines are 
Baraffe et al.l ( 1998h . The bigger symbols with error bars indicate the positions 



of the 9 objects showing water-vapor absorption (squares indicate the 3 objects whose photometry might be 
contaminated by nearby bright /saturated stars). The masses implied for the objects showing water- vapor 
absorption bands are approximately in the range 0.2-1 Mq. 
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served ^0.002 deg'^ area is <1. In addition, 
active galaxies usually show near-IR water 
lines in emission becaus e of the presence o f 
water vapor masers (e.g. Braatz et al.ll2003 ). 
Hence, they are not a likely source of con- 
tamination. 

Artifacts: Another hypothesis to check is 
whether there is some artifact in the data 
that has biased us towards finding cluster- 
centric sources. Our images are not sig- 
nificantly affected by instrumental effects, 
such as vignetting, optical ghosts or image 
degradation moving away from the center, 
that could have caused bias in the selec- 
tion. All 9 selected sources are well de- 
tected in our images, with a signal-to-noise 
ratio better than 5 in all 6 filters (Fig- 
ure [8]). Moreover, all of them are visi- 
ble, though very faint, in the ISAAC and 
ACS images (Sect. |4]). This demonstrates 
that they are not spurious detections. We 
have also checked whether the magnitudes 
of these 9 objects could be somehow affected 
by nearby bright/saturated stars. Also this 
possibility appears to be unlikely for var- 
ious reasons: a) as explained in Sect. [21 
our observing strategy ensures that persis- 
tence spots due to saturated sources do not 
compromise the photometry; b) we excluded 
from our analysis the innermost 5"region of 
the cluster, where there is a high concentra- 
tion of bright OB stars that are saturated; c) 
only 3 (see ID 4, ID 8 and ID 9 in Figure [S]) 
out of the 9 sources are located within a few 
pixels of bright/saturated sources that might 
have contaminated their photometry, while 
all the others are not; these 3 objects are 
marked in Table [S] and in Figure [1] and [5][71 
We adopted PSF fitting magnitudes and lo- 
cal background estimates (see Sect.[2|), which 
should minimize photometric contamination 
even in those 3 cases. However, even in the 
most conservative case in which we exclude 
these 3 objects, we remain with 6 clustered 
objects showing water absorption, a number 
still hard to explain invoking field contami- 
nation, as seen above. 

Binary systems: The 9 objects could be 
BDs belonging to NGC 3603 which are de- 



tected because they are in binary systems 
and, hence, more luminous. Because the 
BD population of NGC 3603 is stih not 
explored, we cannot predict the number 
of binary BDs in the cluster. However, 
they are to be expected, because the frac- 
tion of BDs comprised in binaries mea- 
sured in the field is of the order of 15% 
for a primary mass below 0.1 Mf^ (Lada 



2006HThies fc Kroupa|[2007HBate|[20"09t ) and 



appears to be even higher in star form- 
ing regions and young clu sters of age simi- 
lar to NGC 3603 (~20%; iKononackv elal 



20071 and refer e nces therein) . Moreover , 



Burgasser et al. (l2008l) reported the detec 



tion of a field BD, namely 2MASS J093935482448279, 
presenting a highly inflated radius which can 
be reconciled with BD structure models of 
the same effective temperature (600-700 K) 
only assuming that this source is an unre- 
solved, equal-mass binary. 

A BD/BD or BD/planet pair would have a 
luminosity ^0.75 mag brighter than a single 
BD at maximum, e.g. in the case of the two 
BDs of the pair having the sa me mass. Ac- 
cording to the isochrones by iBarafFe et al. 
(Il998l . our 9 objects are about 2 mag 
brighter than what is expected for single 
BDs belonging to the cluster. This esti- 
mate is affected by our photometric uncer- 
tainties and also depends on two uncertain 
parameters: the distance to NGC 3603 and 
the interstellar extinction along the line of 
sight. The photometric uncertainties for our 
9 objects are of the order of 0.10 mag (see 
Table [5]). The cluster distance is known with 



an un certainty of about 1 kpc (jMelena et al 
l2008[ ). i.e. ~0.33 mag in distance modu- 
lus, while the uncertainty on visual extinc- 
tion is /S.Av ~1 mag (see Sect. 13. 2[) . i.e. 
AAj wO.28 mag. Altogether, these uncer- 
tainties might account for an uncertainty on 
absolute J magnitude of ^0.7 mag. Adding 
this 0.7 mag to the ^0.75 mag excess lumi- 
nosity produced by a BD/BD pair, binary 
BDs would be at most 1.5 mag brighter than 
single BDs in the cluster. This difference is 
still smaller than the 2 mag difference seen 
in the CMD (Figure [T]) between the position 
of our candidates and the expected position 



19 



of BDs belonging to the cluster. Thus, the 
hypothesis of these objects being BD/BD 
or BD/planet pairs appears unlikely. To 
justify the observed luminosity, we should 
assume that they are star/BD pairs with 
the primary star having a mass between 0.3 
and 0.6 Mq, so that the BD would be re- 
sponsible for the observed water absorption 
while the primary star would be responsi- 
ble for the observed luminosity. However, in 
this scenario, the luminosity of the primary 
star in the F139M and F153M filters would 
suppress the BD contribution to our non- 
resolved photometry and, hence, we would 
not be able to detect the presence of water 
absorption bands. 

5. A new class of objects ("Bloatars"): 

The alternative hypothesis would be that 
these objects belong to a formerly unrec- 
ognized class of young stars with a sur- 
face temperature similar to those of BDs, 
but with a luminosity characteristic of much 
more massive stars 0.2 — 1 Mq). These 
stars appear to be "bloated", with radii 
far in excess of what standard PMS evolu- 
tionary theory would predict for objects of 
the same temperature, and so we will dub 
them "Bloatars" (bloated stars). All the 
hints we have collected on these 9 objects on 
the basis of our data support this hypothe- 
sis: their cool temperature, their higher lu- 
minosity with respect to field and/or clus- 
ter BDs, their clustered spatial distribution, 
their CMD similar to that of PMS stars in 
NCG 3603. 

In the next sections we explore this possible and 
interesting scenario in more detail. 

5.1. How could "Bloatars" be formed? 

The data presented in this paper suggest that 
we have discovered a new class of PMS stars, des- 
ignated as "Bloatars", with a surface tempera- 
ture similar to those of BDs but with a luminosity 
characteristic of much more massive stars (~0.2- 
1 M0). How could such stars be formed? Es- 
sentially, it would require a recent mass accretion 
event - preferably with added angular momentum. 
Two scenarios are possible. First, that the system 
was formed as a very close binary system, and that 



the secondary has overflowed its Roche Lobe. Sec- 
ond, the Bloatar has recently swallowed an orbit- 
ing hot Jupiter planet. In both these scenarios, 
the merging is driven by tidal interaction, and the 
debris from the merger forms a common envelope 
which is optically-thick to the escape of photons 
from the primary star's photosphere, thus form- 
ing - for a brief period - a "false photosphere" 
with radius much larger, and temperature much 
lower, than the primary star. We will now inves- 
tigate the likely parameters of such tidal mergers, 
and attempt to determine which scenario gives the 
most plausible results. 

5.1.1. Formation Rates and Timescales 

Planetary systems forming hot Jupiters are 
thought to form in two stages, an initial ac- 
cretion period lasting a few Myr, with orbital 
migration driven by dynamical interaction with 
the protoplanetary disk, followed by a chaotic 
era lasting as long as 100 Myr in which inter ac- 
tions between the newly formed planets dominate 



(|Juric fc Tremaine Il2008[ ). 



The age of the NGC 3603 clust er stars varies 



from z ero up to more than 10 Myr. iBeccari et al 



(|201fl[ ) show that the star formation rate has been 
increasing rapidly over the past few Myr, with the 
bulk of the stars forming over the past ~2-3 Myr. 
Thus, the age of the cluster stars is comparable 
to the length of the first phase of evolution of the 
planetary system. Any planet which is going to 
merge with its star must have therefore already 
migrated to an orbit in which continued tidal in- 
teraction will cause it to inspiral in a timescale of 
less than about 1-3 Myr. 

The inspiral time is determined by a dimen- 
sionless parameter Q representing the ratio of the 
binding energy and the increase in binding energy 
over one orbital period. This is driven primarily 
by tidal dissipation in the satellite body, but may 
also be augmente d by magnetic braking a nd spin 



orbit alignment ([Barker fc Ogilvie 



cally, 10^ < Q < 10*^ (iPeale 19991). The inspi- 
ral timescale is given by 1 
iLevrard et"aLll2009[) : 



2009h . Typi- 



Dobbs-DbMn et al.ll2004t 



2Q (M, 



Win V ni 



(4) 



where n is the orbital angular frequency, the 
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Fig. 8. — Images of the 9 objects showing water-vapor absorption in each on the WFC3/IR fiher used in 
this work, as indicated in the labels. Each snapshot covers an area of 5.2"x5.2"; north is up and east to the 
left. 21 



star has mass and radius i?, , and the planet 
has mass rUp and is orbiting with semi-major axis 
a. If we want the star to inspiral in ~ 1 Myr, 
then the hot Jupiter must initiahy be at a ra- 
dius of order 0.02 - 0.03 AU (3 - 4 x 10"cm) 
with an orbital period of the order of a day. Hot 
Jupiters with such short orbital timescales have 
indeed been observed, e.g. t he eclipsing system 



WASP-18 (,Hellier et al.ll2009l ). If the companion 



was not a planet, but a star in its own right, the 
infall timescale from this initial radius would be 
perhaps a factor of 100 times longer. 

The number, Nb, of Bloatars we expect to see is 
Nb ~ fNc^ where Nq is the number of stars in 
the cluster with mean age tc ^ Myr, / is the 
fraction of stars having hot Jupiters which can 
inspiral in the lifetime of the cluster, and Ib is 
the lifetime of the bloated phase, before the ac- 
creted envelope once again becomes optically thin 
to the escape of photons from the true stellar 
photosphere. According to the most recent esti- 
mates, though still uncertain, the value of / is 
expected to be within the range 0.1-0.2 (see, i.e., 
[Cumming et al . .200 J), while Nr. in NGC 3603 is of 
the order of ~ lOOO ^jfieccari et al.l20lA and Nb is 
^ 10 according to this work. Thus, ^ 5 x 10* yr. 
Bloatars must be short-lived, and therefore com- 
paratively rare objects. 

5.1.2. Orbital parameters & Surface tempera- 
tures 

An inspiralling planet or companion star will 
start to be tidally disrupted w hen it comes to fi ll 
its Roche lobe. According to iPaczvnskil (|l97lf ). 

for mass ratios in the range < < 0.8 this 

will occur when the radius of the companion, r-p, 
exceeds a critical radius, i?cr, given by: 



> Rc. 



0.46224a 



1/3 



(5) 



We assume that the tidally stripped debris from 
the companion form a fast-rotating common enve- 
lope with a radius i?2 — {a + r-p). The stripped 
debris form a "false photosphere" to the star; if 
the tidal debris are optically thick, the effective 
radius of the parent star is '~ R2 and therefore the 
effective temperature of the "false photosphere" is 
given by: 



off 



1/2 



(6) 



Presumably, the effective temperature would 
be higher if viewed from the polar direction, be- 
cause the common envelope would be ellipsoidal 
in shape. If the companion is a lower mass star, 
then nuclear burning is not likely to be immedi- 
ately extinguished, and the effective temperature 
is somewhat higher than that given by equation [B] 
owing to the additional luminosity of the compan- 
ion star circulating within the common envelope. 

Equations [5] and [5] provide the basis for deter- 
mining the most likely formation scenario for these 
Bloatars. Table [7] summarizes the computed effec- 
tive temperature of the bloated photosphere for 
different combinations of mass/temperature of the 
primary star and mass/orbital radius of the com- 
panion planet and shows that a companion hot 
Jupiter provides a much more likely match to the 
observations. 

Indeed, if the companion is a star in its own 
right, the break-up radius is smaller than if it was 
a hot Jupiter. For a companion of mass of 0.21 
orbiting around a primary star of 0.4 — 1.0 Mq, 
we get 1.35 < R2/RQ < 1.56 and 2800 < T^s < 
4700K. At the other extreme, for a companion of 
mass of 0.66 orbiting around a primary star 
of 0.8 - 1.0 Mq, we get 2.8 < i?2/i?0 < 2.9 and 
3000 < Toff < 3600K. These temperatures are too 
high to explain the observations. This fact cou- 
pled with the long inspiral timescale of such mas- 
sive companions makes it most unlikely that the 
Bloatars are simply common-envelope binaries. 

Even a planet like WASP-18 ( nip ~ lOMj - 
lO'^M©; Tp - i?j - O.14i?0) would likewise pro- 
vide too high an effective temperature in the com- 
mon envelope. In the stellar mass range 0.5 — 1.0 
Mq the computed common envelope effective tem- 
perature lies in the range 2800 < Tcs < 4700K. 

For a standard "hot Jupiter" ( nip ^ 1.0 Afj ~ 
lO^^M©; r-p - i?j - O.14i?0) the computed 
common envelope effective temperature lies in the 
range 1970 < Toff < 3300K over the stellar mass 
range 0.5 — 1.0 Mq. However, the tidal heating of 
the companion planet is likely to bloat the planet 
up to larger radius, causing it to fill its Roche Lobe 
earlier and at greater radius from the primary star 
( Leconte et al. 2010f ). If the radius of the planet 
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Table 7: Expected effective temperature of the bloated photosphere (Te/y), according to Equation [51 for 
different combinations of mass/temperature of the primary star (M^,, T^) and mass/orbital radius (m^, r^) 
of the companion planet. 





mp (Mq) 
rp (Ro) 


0.001 
0.14 


0.01 
0.14 


0.001 
0.25 


0.0005 
0.14 


(Mq) 


(K) 


Te// (K) 










1.00 


5920 




3326 


4755 


2489 


2977 


0.93 


5610 




3075 


4393 


2301 


2752 


0.78 


5150 




2775 


3958 


2076 


2485 


0.69 


4640 




2378 


3388 


1780 


2130 


0.47 


3920 




1970 


2794 


1474 


1765 


0.21 


3120 




1266 


1775 


948 


1137 



is increased to Vp ^ O.2i?0, then the computed 
common envelope effective temperature lies in the 
range 1470 < Tcff < 2490K over the stellar mass 
range 0.5 — 1.0 M0 . This nicely brackets the range 
of parameters we derive for our nine stars (0.5 — 1.0 
Mq and 1700 < T^s < 2200K). 

Finally we ask: is the envelope formed suffi- 
ciently dense to form a "false photosphere" to the 
central star? If we adopt the bloated hot Jupiter 
parameters of the previous paragraph, then the 
common envelope will have a local density in the 
range 3 — 6 x 10^^ gm cm~^ and a column density 
~ 10^ gm cm^^. Thus, although the atmosphere 
formed is tenuous, it has more than enough of a 
column density to be extremely optically thick to 
the escape of radiation. 

5.2. Why have "Bloatars" escaped detec- 
tion so far? 

The idea of stars "swallowing p lanets" is not 
new in the literature. For example, Siess &: Livid 
(|l999al ) presented simulations of the accretion of 
massive planets or br own dwarfs by asyni ptotic 
giant branch stars. In Siess fc Livid (|l999bl) they 
also extended this concept to solar-mass stars 
on the red giant branch. Moreover, a few gi- 
ant and dwarf stars have been suggested to un- 
dergo atmospheric expansion and swallowing of 
relat ively massive plane t s in close orbits (see, 
e.g., Israe lian et al. 2001: 'Rett er fc Maronj 12003 : 
Li et al.ll2010; .Fossati et al..,20l'o[l 

On the other hand, all these detections involve 
stellar systems in later evolutionary phases than 
our 9 Bloatar candidates, which are only a few 



Myr old. Thus, a very important issue we need 
to address is why Bloatars have not been detected 
so far in the many photometric and spectroscopic 
surveys of young clusters that are much closer than 
NGC 3603, i.e. in the solar neighborhood (d<l-2 
kpc). 

The probability of detecting a Bloatar is ex- 
pected to be very small. As we saw in Sect. 15. ll the 
"bloated" phase would occur at ages comparable 
with the first phase of evolution of the planetary 
system (1-3 Myrs) and is expected to last only for 
about 5 X 10* yr. In NGC 3603, whose rich PMS 
population and youth provide the circumstances 
most favorable for the detection of Bloatars, they 
would represent just 1% of the ^1 Myrs old PMS 
population. 

These considerations imply that Bloatars could 
have been found only in nearby star forming re- 
gions or clusters younger than a few Myrs. The 
probability of detecting a Bloatar in this restricted 
sample is much smaller than in a massive cluster 
like NGC 3603 for several reasons: 

• Nearby 1-3 Myrs old star forming regions 
and clusters have a typical number of con- 
firmed P MS sta rs of the order of 10^ rather 
than 10^ (jReipu rth 2 008a.b. ). Thus, if we as- 



sume that Bloatars represent just 1% of the 
PMS population, only a few of them (<4) 
are expected to be present in this kind of 
regions. 

• The spatial extent of nearby star form- 
ing regions, like Orion, Serpens, Lupus is 
of the order of a few square degrees or 
more (jReipurthl l2008allb[ ) and even more 
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compact young open clusters, like IC 348, 
typically extend for more than 0.5 deg^ 



(jMuench et al.l 120071 ) . Thus, spectroscopic 



surveys of the PMS population in these 
regions are often spatially incomplete and 
might have easily missed the few Bloatars 
expected to be present. 

• Photometric surveys for PMS stars and BDs 
are nowadays spatially complete. However, 
as we have seen in Sect. HI broad-band col- 
ors of Bloatars are similar to those of 0.2- 
1 Mq PMS stars. Thus, Bloatars might 
have been confused with more ordinary PMS 
stars in such surveys. Only a handful of re- 
gions have been surveyed combining the use 
of broad and narrow-band p hotometry in 



H2O or CH4 
2000; 'Gorlova' ^200 
Ande rs en ct al.. .2006i) . 



li nes ( s ee, e.g.. iNaiita et al 



Mainzer et al 



2004t 



6. Conclusions 

We have presented the results of a pioneer- 
ing narrow-band WFC3/IR survey of NGC 3603 
with the aim of discovering objects displaying deep 
water-vapor absorption bands, characteristic of 
Brown Dwarfs. The main results of our study are 
as follows: 

• We have developed a method to identify 
BD candidates on the basis of WFC3/IR 
narrow-band imaging in water molecular 
bands. This photometric method provides 
effective temperatures for BDs to an accu- 
racy of ±350K relative to spectroscopic tech- 
niques. This accuracy is not significantly 
affected by either stellar surface gravity or 
uncertainties in the interstellar extinction; 

• The comparison of this calibration rela- 
tion with expectations from models of 
stellar/sub-stellar photospheres supports a 
scenario where both dust gravitational set- 
tling and scattering/absorption operate in 
the atmosphere of very cool objects; 

• Using this technique, we identify 9 objects 
having the colors of BDs, strongly clustered 
towards the luminous core of NGC 3603 but 
far too luminous to be normal BD members 
of this cluster. We argue that it is unlikely 



that these objects are either artifacts of our 
dataset, normal field BDs/M-type giants or 
extra-galactic contaminants and, hence, we 
might have discovered a new class of bloated 
objects having the effective temperatures of 
BDs ("Bloatars"); 

• We explore the interesting scenario in which 
these Bloatars would be young stars that 
have recently tidally ingested a Hot Jupiter, 
the remnants of which is providing a short- 
lived extended photosphere to the central 
star. We predict that, if this scenario is cor- 
rect, these stars would show the signature of 
fast rotation, with the equatorial part of the 
"Bloatar" rotating at the orbital velocity of 
the absorbed planet. 
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